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ABSTRACT 

We present a metallicity analysis of 83 late-type giants within the central 1 pc of the Milky Way. 

K-band spectroscopy of these stars were obtained with the medium-spectral resolution integral-field 
spectrograph NIFS on Gemini North using laser-guide star adaptive optics. Using spectral template 
fitting with the MARCS synthetic spectral grid, we find that there is large variation in metallicity, 
with stars ranging from [M/H] < -1.0 to above solar metallicity. About 6% of the stars have [M/H] 

< -0.5. This result is in contrast to previous observations, with smaller samples, that show stars 
at the Galactic center have approximately solar metallicity with only small variations. Our current 
measurement uncertainties are dominated by systematics in the model, especially at [M/H] > 0, 
where there are stellar lines not represented in the model. However, the conclusion that there are 
low metallicity stars, as well as large variations in metallicity is robust. The metallicity may be an 
indicator of the origin of these stars. The low-metallicity population is consistent with that of globular 
clusters in the Milky Way, but their small fraction likely means that globular cluster infall is not the 
dominant mechanism for forming the Milky Way nuclear star cluster. The majority of stars are at 
or above solar metallicity, which suggests they were formed closer to the Galactic center or from the 
disk. In addition, our results indicate that it will be important for star formation history analyses 
using red giants at the Galactic center to consider the effect of varying metallicity. 

Subject headings: Galaxy: center — stars: late-type — stars: abundances — techniques: high angular 
resolution — techniques: spectroscopic 


1. INTRODUCTION 

The metallicity of stars and stellar populations is an 
important property that allows us to understand their 
formation and subsequent evolution. Metallicity can also 
serve as a signature for separating multiple populations 
of stars formed at different times. When averaged among 
many stars, metallicity can be used to trace star forma¬ 
tion within and between galaxies. 

Chemical abundance measurements of stars in the 
Milky Way have shown that there is a strong gradient in 


metallicity (Feltzing & Chiba 


2013). The metallicity in¬ 


creases from below solar metallicity in the outskirts of the 
Milky Way disk to a bove sol ar metallicity within the cen¬ 
tral 5 kpc (|Lepine et al.|20ll]). While the sample of stars 
in the Milky Way with abundance measurements has 
increased dra matically with spectroscopic surveys such 
as APOGEE ( Nidever et al.| |2014), there are abundance 
measurements of on ly about a dozen stars in the central 
10 pc of the Galaxy flCarr et al. 2000; Ramirez et al. 12000 

lonn^ T)_j tv om /it — Ttt— ——— 


Cunha et al.|2007 |Ryde V Schultheis 2014). These mea- 
surements are consistent with the Galactic t 


trend, with a 
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me an [M /H] = 0T4 =b 0.06, and a dispersion of 0.16 dex 
([Cunha et al.||2007). 

The metallicity measurements of stars within the cen¬ 
tral 10 pc of the Galaxy is important, because they form 
the basis of our interpretation of the formation and prop¬ 
erties of the Milk Way nuclear star cluster. This cluster 
is the most massive (10 7 Mq) in the Galaxy and provides 
us with a template for understanding the nuclei of othe r 
galaxies (Schodel et al.||2009| |Chatzopoulos et al.||2015|) . 
The metallicity measurement helps place the cluster m 
context with the rest of the Galaxy, and serves as a start¬ 
ing assumption when inferring the star formation history 
and initial mass function (IMF) from the infrared lumi¬ 


n osity function (Maness et al. 2007 Pfuhl et al. 2011 
L u et al.| 2013). It is therefore import ant to obtain larger 
spectroscopic samples of stars in this region to obtain ro¬ 
bust measurements their physical properties like [M/H]. 

In this study, we combine high angular resolution spec¬ 
troscopy of red giants with spectral template fitting to di¬ 
rectly constrain their effective temperature, gravity, and 
[M/H]. With adaptive optics (AO) spectroscopy, we are 
able to increase both the depth and the number of stars 
with measured [M/H] in this region by a factor of 8. 
This increase in sample size has revealed a number of 
stars with significantly lower metallicity than have been 
previously measured in this region. We also discuss the 
implications of these measurements. 

2. OBSERVATIONS 


The spectra in this paper were obtained using the Gem¬ 
ini North Near-Infrared Facility Spectrograph (NIFS) 
with the natural and laser guide star adaptive optics sys¬ 
tem ALTAIR. The spectra were obtained using the K 
broad-band filter (1.99 - 2.40 @m; GN-2012A-Q-41 and 
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GN-2014A-Q-71, PI: Do). The observations span a pro¬ 
jected radius of 8 to 22 arcseconds (0.3 to 0.9 pc) from 
Sgr A*. More details about the observations and data 
reduction were presented by (Stpstad et al.||2015|). We 
restrict our analysis here to stars with signal-to-noise ra¬ 
tios (SNR) greater than 35 in order to better utilize weak 
spectral lines, and obtain results that are less sensitive 
to priors on the parameters. We also consider only late- 
type stars (F-type or later) with temperatures between 
2500-7000 K. In total, we analyze 83 stars. 


3. SPECTRAL TEMPLATE FITTING 

We fit the observed spectra to a MARCS grid of syn¬ 
thetic mo dels to obtain physical parameters (Gustafsson 


et al. 2008). The MARCS spectral grid spans a range 
of effective temperature (T e //) between 2500 to 7000 
K, surface gravity (logg) between -0.5 and 5.0 dex, and 
scaled solar metallicity (henceforth described by [M/H]) 
between -4.0 and 1.0 dex. The MARCS grid avail¬ 
able onlin^Jare sampled at intervals of 100 K in T e ff. 
[M/H] is sampled variable steps from 0.25 dex between 
— 1 < [M/H] < 1.0, 0.5 dex between —3 < [M/H] < — 1. 
log# is sampled in steps o f 0.5 dex. We conside r models 
of solar composition from Grevesse et ah] (2007| for this 
analysis. In order to obtain spectra from intermediate 
grid parameters, we use a linear interpolation between 
spectra of neighboring grid points. The MARCS grid 
(R = 20000) is convolved to R = 5400 for comparison 
with the NIFS spectra. The spectral resolution is deter¬ 
mined using the average of the full-width half-maximum 
of isolated OH sky-lines. The spectral resolution vary by 
about 10% between different sky lines though not sys¬ 
tematically with wavelength. This amount of variation 
does not significantly affect the fitted parameters. We 
ignore stellar rotation in our fit because red giants are 
obser ved to have rotational velocities below 10 km s -1 
(|Gray||1989), which is not resolvable at the spectral res¬ 
olution of NIFS. We also limit our fitting range to 2.1 to 
2.291 fi m to excl ude CO lines, which tend to bias our fits 
(see Section 4.2). 


We utilize the Bayesian sa mpler MultiNest (Feroz et al. 
[2009 Buchner et al. 2014) to fit the observed spectra. 
The fit is done by computing the posterior: 


P(6\D) = 


P{D\9)P{9) 

P(D) 


( 1 ) 


where D is the observed spectrum, and the model param¬ 
eters 0 = ( T e ff , log#, [M/H], t^), where v z is the radial 
velocity. The priors on the model parameters are P{0) 
and P(D) is the evidence, which acts as the normaliza¬ 
tion. The combined likelihood for an observed spectrum 
is: 


P(D\9)= II -^exp (-(F x ,obs-F x (0)) 2 /2ei obs ), 

A=A 0 e A,o6 S v27T 

( 2 ) 

where is the observed spectrum, F\(6) is the model 

spectrum evaluated with a given set of model parameters, 
and e\, 0 b s is the 1 a uncertainty for each observed flux 
point. This likelihood assumes that the uncertainty for 

9 http://marcs.astro. mi.se/index.php 


each flux point is approximately Gaussian. For compu¬ 
tational efficiency, we use the log-likelihood in place of 
the likelihood: 

-i A n 

In P(D\9) oc - F ^ 0 )) 2 / e CJ- (3) 

* A=A 0 


We choose to use the MultiNest sampler because the pa¬ 
rameter space is often multi-modal, and we find that 
other techniques like Markov-Chain Monte Carlo are less 
efficient and often converge at a local maximum instead 
of the global best-fit solution. 

In addition, for each model evaluation, we apply a 
least-squares fit for a fourth order polynomial that min¬ 
imizes the difference the observed and model spectrum. 
In this way, the stellar absorption lines are accounted for 
in matching the continuum. This procedure however re¬ 
moves our ability to constrain the effective temperature 
using the shape of the spectrum. 

The priors for the model parameters are chosen based 
on the constraints of the MARCS grid and from stel¬ 
lar evolutionary models for stars that could exist at the 
Galactic center. The prior on T e ff is uniform from 2500 
to 7000 K, appropriate for red giants. The prior on [M/H] 
is uniform from -4.0 to 1.0 dex. To determi ne t he limits 
on log#, we use the PARSEC isochrones (Bressan et al.| 


2012) for ages between 10 6 to 10 10 yrs. We then com¬ 


pared the range of log g with that of the K-band luminos¬ 
ity expected at the Galactic center for these isochrones 
(Figure |T]) . Based on these models, we set the priors 
on log#: -0.5 < log# < 4.0 for K < 12 mag and 2.0 
< log# < 4.5 for K > 12 mag. 


4. CHARACTERIZING UNCERTAINTIES 
4.1. Uncertainty from Interpolation 

While interpolation allows us to produce spectra with 
arbitrary model parameters, it also increases the uncer¬ 
tainty in the model spectrum. In order to character¬ 
ize this uncertainty, we remove a grid point, compute 
an interpolated spectrum at that same point, and then 
fit the stellar parameters of the interpolated spectrum 
using the original grid. This represents the maximum 
deviation due to interpolation, because the interpolated 
spectra used for our analyses will be, at most, about 
half the distance from a reference grid point. Repeating 
this process for the entire grid, we find that cr Teff = 50 
K, g\m/h] = 0.1, and criog# = 0.1, with no systematic 
offsets in the fitted parameters. We therefore include 
the interpolation uncertainty by adding these values in 
quadrature with the statistical uncertainty for the stellar 
parameters. These uncertainties are larger than the sta¬ 
tistical uncertainties, but are small compared to model 
uncertainties. 


4.2. Fit comparison with standard spectra 

In order to assess possible model uncertainties in the 
fits, we derive the physical parameters of stars from the 
SPEX stellar s pectral lib rary for comparison to previous 
measurements (Rayner et al. 2009). While this library 
has lower spectral resolution (R = 2000) than NIFS, it 
is the most complete publicly available spectral library 
in the K-band spanning the range of parameter space 
of our sample. Importantly, many of the SPEX spectra 
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Fig. 1.— Left: PARSEC isochrones (Bressan et al.||2012]) sampled over the range of observed stellar ages at the Galactic center from 
the young stars (< 10 Myr) to red giants (1-10 Gyr). Right: The surface gravity of stars as a function of K magnitude. We utilize these 
relationships to limit the range of logg in the spectral fits. Within the brightness range of most of our measurements (K = 12 - 16 mag), 
surface gravity varies from logg = 4.5 to logg = 2. 


have been observed previously and have ha d their stella r 
parameters measured (tabulated by Cesetti et al. [2013). 
By comparing to these previous values, we can estimate 
the range of model uncertainties. 

We include constraints on log g with knowledge of the 
luminosity of the stars, similar to our analysis of the 
Galactic center stars. Stars of luminosity class III were 
limited to 2 < log# < 4. Stars with luminosity class I 
(supergiants) and II were limited to —0.5 < log# < 2. 
Luminosity class V (main sequence) were limited to 
3 < log# <5.5. We also limit the wavelength range of 
the fit from 2.1 to 2.291 /im, as we find there are signifi¬ 
cant biases in log# and [M/H] whe n including CO lines . 
When compared to the values from Cesetti et al. (2013), 
the mean and standard deviation of the fit residuals are: 
A[m/h] = —0.2, cf[m/h] = 0-3, Aiog ^ = 1.0, cri og ^ = 0.9, 
A T eff = 50 K, and OT e * f — 400. Figure [3] shows the 
correlation between our fits to that of the literature, and 
some examples of SPEX spectra compared to our sample 
with similar [M/H]. The offsets from the reference mea¬ 
surements are likely due to systematics in the MARCS 
model spectra, and the fact that the previous measure¬ 
ments were made by different authors, which can differ 
often by 0.2 dex and, in some cases, up to 1 dex in [M/H] 
dCesetti et al.|2013j). Observation of these spectral stan¬ 
dards" with NIF'S will be useful to further quantify the 
comparison between different abundance measurement 
techniques. We include these rms values in the uncer¬ 
tainty for our measured parameters by adding them in 
quadrature with the other uncertainties. These system¬ 
atic uncertainties dominate over all other sources. The 
comparison with previous measurements also shows that 
the log# values are prone to biases, but the values for 
[M/H] and T e ff are consistent to the level of 0.3 dex and 
400 K, respectively, which are accurate enough to obtain 
significant constraints on [M/H] and T eff of stars at the 
Galactic center. 



Reference log g (dex) 



Fig. 2. — Comparison of the best-fit val ues for the SPEX spe ctral 
library to literature values tabulated in Cesetti et al.| (2013} for: 
(top) T e ff , (middle) logg, (bottom) [M/A]. More stars from SPEX 
have T e ff measurements than logg or [M/H]. 

4.3. Comparison to the solar spectrum 

We also test our model and fitting method with the 
spectrum of the sun and find a consistent fit (Figure [4]). 
The solar spectrum^] was convolved to a spectral res- 

10 Obtained from http://www.eso.org/sci/facilities/ 
paranal/decommissioned/isaac/tools/spectroscopic_ 
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Fig. 3.— Examples of stars in our sample convolved to the spectral resolution of R = 2000 (black) to compare with a spectrum from the 
bFEX library with a similar [lVl/ fi[ (blue). The best lit MAKUb model (blue) is able to lit most spectral features, borne features, such as 
Mg 1 near 2.10 jam and (Ja 1 near 2.27 u m are not well reproduced by the model. In general, the tits using t he M AKUb grid are consistent 
with previous measurements of [M/HJ for the bPEX library with a dispersion of about 0.3 dex (see Section |T2} 


olution of R = 5400 and fit with the same wavelength 
range as used to fit the NIFS and SPEX spectra ( 2.1 to 
2.291 /im), with uniform priors, T e ff = 2500 to 7000 K, 
log g = -0.5 to 5.0, [M/H] = -2.0 to 1.0, v z = -000 to 
600 km s -i . The best fit values are: T e ff = 5872 K, 
log# = 4.9, [M/H] = -0.089, v z = 0.032 km s -i . For 
comparison, t he fiducial values for th e sun used by the 
MARCS grid (Gustafsson et al. 2008]) are: T e ff = bill 
K, log# = 4.44, [M/H] = 0.00, v z m 0 km s -1 . The these 
values are consistent with the fit given the uncertainty 
from interpolation and es tima tes of the systematic uncer¬ 
tainty derived in Section |4.2| While most of absorption 
lines are well matched, there are a few that are not well 
matched that may be the source of the small systematic 
offsets from the fiducial values. These lines include the 
Fe I at 2.22632 /im, 2.22662 /im, and 2.23990 /im, and 
Ca I at 2.26574 /im. However, excluding spectral fea¬ 
tures with the large deviations have negligible impact on 
the fitted parameters, as there are many spectral features 
that are incorporated in the global fit (see Appendix for 
more details). 


5. RESULTS 

For each of the 83 spectra, we fit for the 4 physical 
parameters: T e //, log//, [M/HJ, and v z . We report the 
central value of the probability distribution for each of 
these parameters marginalized over all other model pa¬ 
rameters in Table 1. For most stars in the sample, there 
are correlations between T e /f, log#, and [M/HJ, depend¬ 
ing on the star, which emphasizes the necessity to fit 
them simultaneously. As expected, there is no correla¬ 
tion between v z and the other 3 parameters. Figure [5] 
shows an example of the joint probability distribution 
functions (PDFs) for the 4 stellar parameters in the fit 


standards.html using NSO/Kitt Peak FTS data produced 
by INbF/iNOAO. 


for the star NE1-1 003. These PDFs show that the sta¬ 
tistical uncertainties are generally very small compared 
to the systematic uncertainties estimated in Section [4j 
The fits to the physical parameters reveal the existence 
of stars with low metallicity at the Galactic center. These 
stars have unusual spectra compared to the rest of the 
sam ple, showing low CO, Na, and Ca equivalent widths 
(e.g. Pfuhl et al. 2011 jDo et al.||2013 ). In Figure [ 6 j we 
show selected stars with a range of metallicity, from low 
to super-solar. For comparison, we also show the best- 
fit spectra when the metallicity is fixed to [M/H] = 0.0. 
Certain lines such as Si I and Fe I are sensitive to [M/H], 
while others, especially H I at 2.1661 /im, are more sensi¬ 
tive to temperature. For example, the star NE1-1 003 is 
better fit with low metallicity of [M/H] = -1.28 and T e ff 
= 4136 K, compared to one with fixed [M/H] = 0.0, and 
best fit T e ff = 5157 K (Figure [ 6 ). NE 1-1 003 lacks the 
stronger Br 7 line at 2.1661 /iuirhat would be required 
for the T e ff = 5157 K fit, and contains weak Fe I and Si 
I lines that are more consistent with low [M/H]. Some of 
the lines that show mismatches between the model and 
observations, such as Fe I at 2.23990 /im, and the Ca I 
line at 2.26574 /im (Figure [ 6 ] and [7]). These lines are also 
mismatched in our solar spectrum comparison (Section 


4.3). These mismatches may be a concern if deriving 
individual elemental abundances, but as our fit aims to 
determine relative scaled solar abundances, the impact of 
these mis matc h is lesse n, as verified by our comparisons 
in Section 4.2 and 4.3 While there may be systematic 
uncertainties in the absolute measurement of [M/H], the 
conclusion that these stars must be low metallicity com¬ 
pared to most of the sample is robust. In total, 5 of 83 
( 6 %) of the stars have [M/H] < —0.5 (Figure [ 8 ]). 

Most of the stars have higher than solar metallicity, 
with a sample mean of [M/H] = 0.4 and a standard de¬ 
viation of 0.4 dex. This suggests that there are many 
stars with super-solar metallicity at the Galactic cen- 
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Fig. 4.— The best fit MARCS model (blue) compared to the observed solar spectrum in the near-infrared (black; NSO/NOAO). The 
fiducial solar model is also shown (green). In general most of the lines are well fit, with the best fit within the range of uncertainties we 
find from comparison to the SPEX spectral templates. A few lines have are not well matched by the model, which likely is a source of 
systematic uncertainty in our fits. 



Teff(K) logg [M/H] Vz (km/s) 


Fig. 5.— The joint probability distributions for T e ff, logg, [M/H], 
the stellar parameters for NE1-1 003. 

ter, though they are likely subject to greater systematic 
uncertainties. The MARCS models are unable to repro¬ 
duce many features in the metal-rich sample. For exam¬ 
ple, in Figure [6j the star E5-2 001 shows strong features 
at 2.1898 and 2.2653 fi m, which are not present in the 
MARCS models. This may be due to differences from so¬ 
lar composition or missing atomic and molecular lines in 
generating the spectra. In addition, at medium spectral- 
resolution, the difference between solar-metallicity and 


and v z , as well as their 1-dimensional probability distributions for 

high-metallicity spectra are much smaller than that be¬ 
tween solar-metallicity and low-metallicity sources (Fig¬ 
ure^. In order to resolve these issues, it will be neces¬ 
sary to expand the comparison to non-solar composition 
models, and to obtain a larger sample of high spectral 
resolution observations to calibrate these [M/H] measure¬ 
ments. While the large spread in [M/H] is likely robust 
as they are relative measurements, conclusions relying 
on the location of the peak or shape of the metallicity 
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Fig. 6. — Examples of spectra with a range of [M/H] along with their best fits (blue) compared to a fit with [M/H] fixed to 0.0. Labeled 
are the best fit [M/H] values. The three stars show examples of metal-poor, solar metallicity, and super-solar metallicity stars. These two 
wavelength regions of K-band were chosen to illustrate the combination of temperature-sensitive lines such as H I compared to metallicity 
sensitive lines such as Fe I and Si I. The MARCS spectra appear to a good match at solar or below solar metallicities. At high metallicities, 
the model is unable to reproduce some of the observed lines, likely leading to larger systematic uncertainties in the value of [M/H] for these 
stars. 


distribution are probably unreliable at this time. 

6. DISCUSSION 

6.1. Foreground/background sources and spectroscopic 
contamination 

Interpretation of the results of this study depends 
in part on whether the sources belong to the nu¬ 
clear star cluster, or are foreground or background 
sources. Some obvious foreground sources can be ex¬ 
cluded based o n their blu e colors , and are not in the 
present sample (Stpstad et al. 2015). The remaining fore- 
ground/background contaminants are likely stars from 
the inner bulge of the Milky Way. The most recent 
estimate of the number density of stars i n th e _i_nn£r 
bulge in the near-infrared was made by |Clarkson et al7| 
(2012) with their measurement of the proper motion of 
the Arches star cluster, about 26 pc in projection from 
the Galactic center. Using their proper motion identifica¬ 
tion of field stars, we find that the likely number of fore¬ 
ground/background sources is 0.08 stars/arcsec -2 with 
K < 15.5. With a coverage of about 99 arcsec 2 in the 
current survey, we expect about 8 stars to be from the 
inner bulge. It is unlikely that the low-metallicity stars 
found in this study belong to the inner bulge, as the in¬ 
ner bulge was measured to have [M/H] = —0.16 =b 0.12 
(|Rich et al.||2012| |Ryde & Schultheis||2014 ). 

A second source contamination, especially at low- 
metallicity, may be Milky Way halo interlopers. Using 
the Besancon Model (Robin et al. 2003), we estimated 
the expected number oFTTalo stars expected between 10 


- 18 mag in K (including extinction) using an extinction 
law of Ay = 3.5mag/kpc. Within our field of view, the 
model predicts 0.0004 halo stars, which is negligible. 

The presence of a nearby companion would affect the 
spectra of the stars, biasing the measured stellar pa¬ 
rameters. Early-type stellar (B-type or earlier) compan¬ 
ions will introduce the strongest bias. These compan¬ 
ions may either be physically close binaries or projected 
pairs. However, this is unlikely for two reasons: (1) there 
are ver y few 2) ea rly-typ e stars in the region of this 
survey (Stpstad et al. 2015), so chance superposition is 
negligible, (2) physical binaries of a late-type giant and 
an early-type star are unlikely based on stellar evolution, 
as the early-type star will have a lifetime of < 100 Myr, 
while the late-type giant is closer to 1 Gyr in age. 

6.2. Comparison with previous work 

There have been few measurements of individual 
[M/H] values for stars in the central parsec. About a 
dozen bright late-type stars, most of which are red super¬ 
giants or AGB stars ha ve been measured at high spectral 
resolution. Carr et al. (2000) was the first to measure the 
metallicity of the red supergiant, IRS 7 ([M/H] = -0.02 
±0.13), located abou t 5.5" (0.22 pc) from the Galactic 
center. Subsequently, |Ramfrez et al7| (2000) measured 9 
red supergiants and AGB stars (K < 8 mag) within 2.5 
pc of the Galactic center and found a mean [M/H] = 
0.12±0.22. The [M/H] of their sample range from -0.29 
to 0.49, wit h typical uncertainties of 0.3 dex. [Cunha 
et al. (|2007) re-observed 5 of the same sources and found 
a mean metallicity of [M/H] = 0.14, and a dispersion of 
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Fig. 7.— Comparison of the best fit model for NE1-1 003 ([M/H] = -1.28) to [M/H] — -0.5 and 0.0 to show the effect on the spectrum 
with increasing metallicity. 



[M/H] 

Fig. 8. — Distribution of [M/H] measurements for stars with 
[M/H] < 0.3 (remaining stars have higher measured [M/H], but 
are likely subject to larger systematic uncertainties). 


fainter than the previous measurements with the bulk of 
the sample between 12 < K < 15.5 mag. These stars 
are thus closer to the red clump than the more evolved 
red-supergiants and AGB stars previously observed in 
this region. In addition, there may be systematic ef¬ 
fects introduced by our medium-resolution spectra com¬ 
pared to the previous high-spectral resolution observa¬ 
tions. The method of deriving [M/ H] also vary between 
the different studies. For example, dRamfr ez et al.l f2000) 
used the MOOG spectral synt hesis code ([Sneden 1973) 
while Ryde & Schultheis (2014) used Spectroscopy Made 


Easyl^Ve ilenti & Piskunov] 1996). In order to characterize 
these uncertainties, it will be useful to obtain compare 
these methods for stars at the Galactic center. If con¬ 
firmed, some of these stars may be the most metal-rich 
ever found in the Galaxy. 


0.16 dex. More recently, Ryde & Schultheis (2014) ob¬ 
served 9 fainter giants with K < 12 mag within 2m and 
3.5' (6 to 8.4 pc) from the Galactic center. These stars 
have [M/H] = 0.11T0.15, similar to the observations in 
the central parsec. 

The measured range of [M/H] in our current sample 
is large compared to previous measurements. The value 
of [M/H] for our sample of 83 stars range from < —1.0 
to > 0.8 dex. This wide range may be the result of 
the larger sample of stars, which allows us to detect the 
rarer low-metallicity stars. Our sample does not over¬ 
lap the previous ones, and the stars are also significantly 


6.3. Implication for the formation of the nuclear star 

cluster 

The metallicity of the nuclear star cluster encodes in¬ 
formation about its history and initial chemical composi¬ 
tion. Until this study, no low-metallicity stars had been 
found there. They may represent stars that were in this 
region early in the history of the Milky Way, or may 
have migrated over time to the Galactic center. Another 
possibility is that these stars arrived through the infall 
of globular clusters, which we see today with average 
[M/H] similar to the low [M/H] stars detected in this 
study. In fact, one of the prevailing theories for nuclear 





























of globular clusters over time (e.g. 

Tremaine et al. 

1975; 

Capuzzo-Dolcetta & Miocchi 2008 

Antonini et alG 

MM 

Measurement of the metallicity 0 

t these stars offers a 


new opportunity to test this theory. Globular clusters 
in the Milk y Way^ generally have [M/H] of about -1.0 
to 


-2.0 (H 


arris 


1996). The highest metallicity clusters 


are slightly below solar metallicity, such as Terzan 5 at 
[M/H] = -0.23. If the globular clusters we see today are 
representative of the past population, the small fraction 
of low metallicity stars in our sample suggests that glob¬ 
ular clusters are only a small contributor to the origin of 
the Milky Way nuclear star cluster. 

Accounting for the existence of super-solar metallic¬ 
ity stars will be important to determine the origin of 
the nuclear star cluster, as these represent the major¬ 
ity of stars. The source of these stars may be from the 
Galactic disk, where the metallicity of stars are generally 
higher. Super-metal rich stars in the solar neighborhood, 


Galactic bulge (Trevisan et al. 

2011 

Bensby et al. 2013; 

Feltzing & Chiba|2013p. In order to 0 
of these stars in the nuclear star clus 

etermine the origin 

ter, we will need to 


identify the systematic uncertainties in the model (Fig¬ 
ure [6]), and obtain high resolution spectroscopy to mea¬ 
sure stellar abundance ratios. Abundance ratios such as 
[a/Fe], [Ca/Fe], [Si/Fe], etc., can be matched to signa¬ 
tures from other parts of the Galaxy such as the disk, 
bulge, globular cluster, or local dwarf galaxies. 

6.4. Implication for stellar population analyses 

The detection of the large spread of metallicity of stars 
at the Galactic center has the potential to change the way 
we study this region. It may be necessary to revisit the 
measurements of the sta r formation history and IMF of 
the nuclear star clus ter ( Maness et al.||2007 


2011 Lu et al.|2013 ). These studies have 


Pfuhl et al. 


_so far assumed 

solar metallicity in their models. These analyses require 
a translation from the luminosity function into a mass 
function through the use of evolutionary and atmospheric 
models, which can strongly depend on the metallicity of 
the stars. In Figure [9j we show a comparison between the 
iGband luminosity function resulting from 0.5 to 2 M 0 
stars of a cluster with an age of 10 9 yrs with [M/H] = - 
1.0, 0.0, and +0._5jising the PARSEC evolutionary tracks 
(Bressan et al. 2012). These isochrones show that there 
are significant differences in the luminosity and tempera¬ 
ture of stars at different metallicities, which will need to 
be accounted for when deriving a star formation history 
of this region. 

7. CONCLUSION 

We presented stellar parameter estimates using K-band 
spectra of 83 stars within 1 pc of the center of the Galaxy, 
and we find a significant spread in metallicity of stars, 
ranging from 10 times below solar to super-solar metal¬ 
licities (about 6% of the sample have [M/H] below -0.5). 
This variation in metallicity shows that the Milky Way 
nuclear star cluster is not composed of a simple stellar 
population, which previous work has also shown with 
spectroscopy and luminosity fu nctions (e.g. Maness et al. 


2007 Chatzopoulos et al. 2015). Future measurements 
of the star formation history of this region will need to 
take the variation in metallicity into consideration. The 



10000 9000 8000 7000 6000 5000 4000 3000 2000 

Temperature (K) 


Fig. 9. — Evolutionary tracks from the PARSEC model ( Bressan] 
|et al.|[20121) at [M/H] = -1.0, 0, +0.5 at 10 9 yrs, showing the 
difference m temperature and AT-band magnitude. The AT-band 
magnitudes are adjusted for the distance to the Galactic center 
(8 kpc) and the AT-band extinction (Ak = 3.0). This shows that 
interpretation of the stellar population of the Galactic center can 
change significantly based on its metallicity. 

low-metallicity stars found in this study are consistent 
with the range of [M/H] observed for globular clusters, 
which is predicted by the theory that infalling globular 
clust ers contribute to the build up of nuclear star clus¬ 
ters d Capuzzo-Dolcetta fc Miocchi||2008l |Antonini et al. 
2012|) This fraction however is small, and thus this 


mechanism is not likely to be the source of the Milky 
Way nuclear star cluster. The evidence for super-solar 
metallicity, on the other hand, points towards a signifi¬ 
cant contributi on from stars with origi ns near the center 
of the Galaxy (jFeltzing & Chiba 2013). 

This paper uses data from Gemini observations GN- 
2012A-Q-41 and GN-2014A-Q-71. We wish to thank 
the anonymous referee and Norm Murray for useful com¬ 
ments. 
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APPENDIX 

TESTING THE EFFECT OF MISMATCHED SPECTRAL FEATURES 

Here we analyze possible effects on the stellar parameter fits from mismatches of some spectral features from the 
model spectra. For example, in the comparison with the Solar spectrum (Figure [4]), there are a few lines, such as 
the Fe I at 2.22632 /im, 2.22662 /im that are not well represented in the model. To determine the effect of this 
mismatch on the fitted parameters, we compute the residuals between the best fit model and the Solar spectrum to 
determine regions that deviate by more than 3%. 5 features exceed this level of deviation (Figure [l0|. We exclude the 
4 spectral channels around wavelength regions at 2.17688, 2.19012, 2.2386.9, 2.2263.4 /im, and lo spectral channels 
around 2.12418 /im, and refit the spectrum. The best fit parameters are T e ff = 5902 K, log// = 4.9, [M/H] = -0.07, 
v z = 0.45 km s -1 . These values are consistent with the fit to the full spectrum ( T e ff = 5872 K, log// = 4.9, [M/ H] 
= -0.089, v z = 0.032 km s _1 ). We also repeat the analysis of the SPEX spectral library described in Section |4.2| by 
masking the same features. We find that the mean difference and standard deviation between the sample with clipping 
and the one without to be: T e ff = —2.6T26 K, log// = —0.010T0.14 dex, [M/H] = 0.002T0.04 dex. These differences 
are about an order of magnitude smaller than the statistical uncertainties, so we conclude that these spectral features 
do not affect the fitted stellar parameters. 


TABLE 1 

Measured physical parameters of stars at the Galactic center 


Name 3, 

K 

SNR 

RA Offset (") 

DEC Offset (") 

T eff (K) 

° T eff b (K) 

log 9 

<Tl 0 g g h 

[M/H] 

^[M / H ]° 

v z (km s *■) 

E5-1 001 

12.0 

42 

15.17 

-0.49 

3497 

413 

3.04 

0.91 

0.96 

0.32 

-54 

E5-1 002 

12.6 

54 

14.52 

1.47 

3671 

414 

2.84 

0.91 

0.55 

0.32 

-64 

E5-1 003 

13.2 

44 

15.00 

1.06 

3597 

414 

3.09 

0.91 

0.85 

0.32 

-58 

E5-1 006 

14.7 

40 

16.80 

0.36 

4076 

417 

3.86 

0.91 

0.44 

0.32 

-29 

E5-1 007 

15.4 

41 

15.25 

-0.23 

3754 

414 

3.18 

0.91 

0.40 

0.32 

-57 

E5-1 008 

15.1 

36 

15.02 

0.68 

3915 

418 

3.56 

0.91 

0.43 

0.32 

-112 

E5-1 009 

14.8 

39 

13.83 

-0.29 

4511 

419 

3.91 

0.91 

0.40 

0.32 

-61 

E5-1 012 

15.2 

36 

14.12 

-0.04 

4072 

420 

3.64 

0.91 

0.63 

0.32 

112 

E5-1 014 

15.4 

38 

16.53 

1.32 

4197 

421 

3.82 

0.91 

0.27 

0.33 

-52 

E5-1 015 

15.7 

39 

14.78 

-0.36 

4180 

415 

3.91 

0.91 

0.25 

0.32 

-75 

E5-1 016 

15.5 

42 

15.75 

-0.69 

4145 

414 

3.89 

0.91 

0.13 

0.32 

58 
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2.10 2.15 2.20 2.25 

Wavelength (micron) 


Fig. 10. — Residual of the best fit MARCS model to the Solar spectrum. In order to test the robustness of our fit against large deviations 
in the model, we exclude regions where the residuals are greater than 3%. Excluding these regions do not significantly impact the fits 
for the solar parameters ([M/H] changed by less than 0.02 dex, an order of magnitude smaller than the uncertainties for the stars in our 
sample). 

TABLE 1 — Continued 


Name 3, 

K 

SNR 

RA Offset (") 

DEC Offset (") 

Teff (K) 

° T ef f b (K) 

log# 

<Rog g h 

[M/H] 

a [M/H ] b 

v z (km s x ) 

E5-1 026 

16.2 

36 

15.61 

0.30 

4310 

418 

3.93 

0.91 

-0.05 

0.33 

-134 

E5-1 042 

14.9 

35 

15.28 

1.05 

4161 

417 

3.82 

0.91 

0.56 

0.32 

228 

E5-2 001 

11.3 

60 

15.47 

-4.21 

3519 

413 

2.92 

0.91 

0.73 

0.32 

-63 

E5-2 003 

14.0 

49 

13.76 

-2.86 

3932 

419 

3.49 

0.91 

0.02 

0.32 

-189 

E5-2 004 

13.9 

51 

15.01 

-3.50 

3808 

414 

3.29 

0.91 

0.39 

0.32 

-13 

E5-2 005 

14.2 

38 

14.19 

-2.05 

4328 

415 

3.85 

0.91 

0.40 

0.32 

-29 

E5-2 006 

14.3 

36 

13.29 

-3.38 

4071 

418 

3.87 

0.91 

0.45 

0.32 

157 

E5-2 009 

14.9 

38 

14.49 

-2.26 

4251 

423 

3.09 

0.91 

0.00 

0.33 

4 

E5-2 010 

15.1 

51 

15.67 

-2.25 

4434 

416 

3.91 

0.91 

0.18 

0.32 

-37 

E5-2 Oil 

15.4 

41 

15.72 

-3.19 

4613 

415 

3.94 

0.91 

0.16 

0.32 

-13 

E5-2 019 

15.6 

37 

15.10 

-3.87 

4131 

415 

3.69 

0.91 

0.11 

0.32 

-21 

E5-2 020 

15.7 

36 

14.65 

-2.91 

4066 

420 

3.79 

0.91 

0.42 

0.32 

-82 

E6-1 001 

12.0 

53 

17.88 

0.89 

3687 

413 

3.17 

0.91 

0.60 

0.32 

84 

E6-1 002 

12.8 

54 

17.95 

-1.80 

3737 

414 

3.19 

0.91 

0.47 

0.32 

-59 

E6-1 003 

13.3 

48 

18.92 

-1.56 

3656 

414 

3.04 

0.91 

0.65 

0.32 

4 

E6-1 004 

13.3 

75 

18.92 

0.21 

4155 

413 

3.68 

0.91 

-0.12 

0.32 

-61 

E6-1 005 

13.8 

58 

17.91 

-1.02 

3974 

414 

3.47 

0.91 

0.33 

0.32 

-218 

E6-1 006 

13.8 

44 

19.45 

-1.53 

3720 

414 

3.25 

0.91 

0.51 

0.32 

-38 

E6-1 008 

14.8 

38 

18.40 

-1.13 

3986 

416 

3.68 

0.91 

0.43 

0.32 

-4 

E6-1 009 

14.9 

55 

17.56 

-0.47 

4358 

413 

3.86 

0.91 

0.41 

0.32 

-124 

E6-1 010 

15.4 

49 

18.18 

-0.47 

4590 

417 

3.63 

0.91 

-0.04 

0.32 

-105 

E6-2 001 

13.0 

66 

19.06 

-2.98 

3859 

414 

3.19 

0.91 

0.43 

0.32 

46 

E6-2 002 

13.5 

50 

16.21 

-3.57 

3759 

413 

3.26 

0.91 

0.39 

0.32 

-52 

E6-2 003 

14.0 

42 

19.07 

-2.56 

3687 

414 

3.07 

0.91 

0.63 

0.32 

1 

E6-2 004 

14.0 

42 

16.54 

-3.67 

3730 

415 

3.08 

0.91 

0.42 

0.32 

-10 

E6-2 005 

14.8 

36 

17.49 

-2.98 

3836 

414 

3.30 

0.91 

0.64 

0.32 

16 

E6-2 006 

14.8 

41 

17.24 

-3.20 

3898 

419 

3.47 

0.91 

0.53 

0.32 

71 

E6-2 008 

15.2 

45 

18.12 

-2.94 

4063 

415 

3.61 

0.91 

0.24 

0.32 

-71 

E6-2 009 

15.1 

42 

17.80 

-4.55 

4422 

416 

3.86 

0.91 

0.16 

0.32 

-49 

E7-1 001 

10.8 

62 

20.46 

-2.29 

3479 

413 

2.92 

0.91 

0.73 

0.32 

-105 

E7-1 002 

11.7 

56 

21.80 

-0.93 

3662 

413 

2.96 

0.91 

0.56 

0.32 

-70 

E7-1 003 

12.1 

41 

21.40 

-1.99 

3591 

414 

2.90 

0.91 

0.76 

0.32 

35 

E7-1 004 

12.3 

45 

20.28 

-1.08 

3594 

414 

2.95 

0.91 

0.72 

0.32 

-120 

E7-1 005 

13.4 

39 

20.95 

-2.11 

3888 

417 

3.41 

0.91 

0.39 

0.32 

23 

E7-1 006 

13.5 

52 

22.73 

-0.77 

3878 

414 

3.29 

0.91 

0.52 

0.32 

-29 

E7-1 007 

13.3 

64 

20.67 

-2.39 

3479 

413 

2.83 

0.91 

0.64 

0.32 

-118 

E7-1 010 

14.6 

39 

20.47 

-1.19 

3663 

416 

2.78 

0.91 

0.41 

0.32 

-135 

E7-1 022 

14.8 

38 

20.98 

-0.17 

4218 

425 

3.63 

0.91 

0.44 

0.32 

38 

E7-2 001 

11.4 

51 

20.78 

-2.97 

3524 

413 

2.85 

0.91 

0.74 

0.32 

-189 

E7-2 002 

14.3 

50 

20.20 

-4.42 

3801 

414 

3.27 

0.91 

0.39 

0.32 

-204 

E7-2 006 

15.4 

56 

20.49 

-3.73 

4313 

416 

3.76 

0.91 

0.17 

0.32 

-90 

Nl-1 001 

13.2 

44 

3.51 

8.17 

3644 

414 

2.99 

0.91 

0.68 

0.32 

-47 

Nl-1 002 

13.3 

81 

1.40 

8.51 

4198 

413 

3.85 

0.91 

-0.91 

0.32 

-72 

Nl-1 003 

13.9 

46 

4.25 

9.00 

3844 

416 

3.36 

0.91 

0.48 

0.32 

-128 

Nl-1 004 

14.1 

37 

3.89 

8.04 

3988 

417 

3.36 

0.91 

0.36 

0.32 

-3 

Nl-1 005 

13.6 

45 

3.50 

9.63 

3658 

414 

3.18 

0.91 

0.75 

0.32 

-122 

Nl-1 007 

14.4 

36 

3.97 

8.89 

3851 

415 

3.22 

0.91 

0.61 

0.32 

79 

Nl-1 045 

12.1 

84 

3.64 

7.22 

4212 

413 

3.97 

0.91 

-1.14 

0.32 

170 

Nl-2 001 

13.3 

56 

5.92 

6.93 

3929 

415 

3.52 

0.91 

0.40 

0.32 

57 

Nl-2 002 

13.3 

51 

6.67 

7.99 

3706 

413 

3.19 

0.91 

0.67 

0.32 

53 

Nl-2 003 

13.5 

41 

5.48 

8.60 

3689 

413 

3.21 

0.91 

0.65 

0.32 

60 



11 


TABLE 1 — Continued 


Name 3, 

K 

SNR 

RA Offset (") 

DEC Offset (") 

Teff (K) 

*T eff h (K) 

log 9 

'Rog g° 

[M/H] 

&[ M / H ] b 

v z (km s x ) 

Nl-2 004 

13.9 

47 

6.33 

7.68 

3844 

413 

3.74 

0.91 

0.36 

0.32 

137 

Nl-2 006 

14.9 

36 

6.08 

8.20 

4204 

420 

3.86 

0.91 

0.39 

0.32 

30 

Nl-2 016 

16.7 

35 

4.86 

7.50 

4484 

416 

3.10 

0.91 

0.13 

0.32 

-83 

N2-1 001 

12.1 

49 

5.79 

11.57 

3593 

413 

3.12 

0.91 

0.69 

0.32 

61 

N2-1 002 

12.1 

110 

4.40 

11.00 

4358 

413 

3.92 

0.91 

-1.06 

0.32 

223 

N2-1 003 

12.8 

43 

3.59 

11.53 

4260 

415 

3.82 

0.91 

-1.20 

0.32 

31 

N2-1 004 

13.1 

37 

3.63 

10.68 

3909 

415 

3.41 

0.91 

0.18 

0.32 

48 

NE1-1 001 

10.4 

54 

11.54 

3.73 

3558 

414 

3.00 

0.91 

0.80 

0.32 

55 

NE1-1 002 

10.7 

52 

8.58 

4.10 

3447 

414 

2.78 

0.91 

0.90 

0.32 

-141 

NE1-1 003 

11.4 

133 

8.61 

3.76 

4125 

413 

3.98 

0.91 

-1.27 

0.32 

-119 

NE1-1 005 

12.4 

50 

9.43 

4.95 

3517 

413 

2.93 

0.91 

0.89 

0.32 

-135 

NE1-1 007 

13.5 

39 

9.68 

2.99 

3710 

414 

3.21 

0.91 

0.65 

0.32 

59 

NE1-1 008 

13.4 

48 

10.54 

4.92 

3625 

415 

3.02 

0.91 

0.81 

0.32 

-106 

NE1-1 009 

13.5 

52 

8.98 

4.25 

3683 

413 

3.27 

0.91 

0.56 

0.32 

-43 

NE1-1 010 

13.9 

43 

10.80 

3.67 

3769 

417 

3.30 

0.91 

0.62 

0.32 

-186 

NE1-1 Oil 

14.2 

46 

9.36 

3.15 

3829 

414 

3.35 

0.91 

0.56 

0.32 

34 

NE1-1 012 

14.5 

41 

11.11 

4.41 

4239 

418 

3.75 

0.91 

0.48 

0.32 

-171 

NE1-1 013 

14.8 

41 

11.08 

4.79 

3885 

422 

3.35 

0.91 

0.33 

0.32 

-225 

NE1-1 014 

15.3 

38 

11.51 

2.86 

4241 

417 

3.92 

0.91 

0.03 

0.32 

54 

NE1-1 018 

15.1 

39 

10.37 

3.57 

4354 

417 

3.79 

0.91 

0.20 

0.33 

-41 

NE1-1 025 

15.3 

52 

9.18 

2.84 

4224 

419 

3.66 

0.91 

0.55 

0.32 

-39 


a Name, K magnitude, and RA & DEC offset from Sgr A* from jStpstad et al.|2015| . 

k Uncertainties include statistical uncertainties, interpolation uncertainties and systematic uncertainties added in quadrature. 








